A possible generalized form of Jarzynski equality 
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The crucial condition in the derivation of the Jarzynski equality ( JE) from the fluctuation theorem 
is that the time integral of the phase space contraction factor can be exactly expressed as the entropy 
production resulting from the heat absorbed by the system from the thermal bath. For the system 
violating this condition, a more general form of JE may exist. This existence is verified by three 
Gedanken experiments and numerical simulations, and may be confirmed by the real experiment in 
the nanoscale. 
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I. INTRODUCTION 

Consider a classic system in contact with a thermal 
bath at constant temperature, and at some time inter- 
val, the system is driven out of the equilibrium by an 
external field. Two groups of equalities are proved to 
still hold for this system. One is the fluctuation theorem 
(FT) 0, S H i, ll which reflects the probability of vi- 
olating the Second Law of Thermodynamics in the non- 
equilibrium process. Another is the Jarzynski equality 
(JE) @ which ensures us to extract the free en- 
ergy difference between two equilibrium states from the 
non-equilibrium work performed on the system in the 
process between these two states. The quantum versions 
P. ITollTlj | and experimental verifications 0, of JE are 
also presented. After it was proposed in 1997, the JE has 
aroused some controversy [Tl 111 111 [rl ITS iTfl El l2l| . 
in which two typical gedanken experiments are quite in- 
teresting and we summarize them as follows. 

Experiment 1: As shown in Fig. imagine that a 
closed container, in contact with a thermal bath at con- 
stant temperature, is divided into two compartments by a 
perfectly thin, frictionless but heavy enough piston, and 
imagine that one compartment initially contains ideal gas 
of N (large enough) particles in equilibrium at tempera- 
ture T, while another compartment is empty. At time t\, 
we remove the pins P\ and Pi , and give the piston a large 
initial velocity v p . The gas will fill the whole container 
with the movement of the piston. After a long time re- 
laxation, the system arrives at an equilibrium state at 
time ti- 

Experiment 2: As shown in Fig. ^3, imagine that a 
closed container, in contact with a thermal bath at con- 
stant temperature, is divided into two compartments by 
a perfectly thin and frictionless plate, and imagine that 
one compartment initially contains ideal gas of N parti- 
cles in equilibrium at temperature T, while another com- 
partment is empty. At time t\, we pull up the plate, and 
the gas will expand and fill the whole container. After a 



long time relaxation, the system reaches an equilibrium 
state at time Here t\ and t% do not require to have 
the same values as those in experiment 1. 
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FIG. 1: Gedanken experiments. 



Assume the initial volume of the gas to be Vi, the 
whole volume of the container to be V% ■ There are two 
common points in the above two experiments: (i) The 
macroscopic work in the expansion process is vanish- 
ing (i.e. W = 0); (ii) After the systems arriving at 
the final equilibrium states, the free energy difference 
is AF = -NTla(V 2 /Vi). The important difference be- 
tween them is that the microscopic work, w, in the first 
experiment is non-vanishing although the macroscopic 
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(w) 



for Vr, 



while w = for 



the second one. Due to this difference, the JE holds in 
the first experiment (i.e. {e~ 0w ) = e~ pAF = V 2 N /V{ 



N 



with /3 = 1/T) |lSl but fails in the second one ((e^™) = 
1, e~@ AF = V^ 1 /Vl 1 ). Jarzynski and Crooks argued that 
the JE fails because the initial distribution function is 
not canonical in the second case [rH |22| . We would like 
to consider this problem from another point of view: the 
initial distribution function is still canonical but a more 
underlying reason makes the JE fail. In other words, 
there is a more general form of JE. The rest of this paper 
is focus on this topic and organized as follows: In Sec.lTTl 
we sketch the derivation of JE from the FT and empha- 
size that the condition of adiabatic incompressibility p3[ 
is crucial to this derivation. In fact, Jarzynski's original 
proof H E3 also requires this condition. In Sec. IIIII we 
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check whether this condition holds or not in the above 
two Gedanken experiments and put forward a generalized 
JE, Ea. (|15(l . The third Gedanken experiment interme- 
diating between the above two experiments is proposed 
and the corresponding numerical simulation verifies the 
existence of generalized JE. In Sec. IIVI we give further 
discussions and a brief summary. 



II. DERIVATION OF JE FROM FT 

An important relation between the FT and JE is that 
the JE can be derived from the FT for time reversible 
stochastic or deterministic dynamics |U H^- Here we 
look through the main idea of Evans' derivation [2^ . 
The phase space of the N particle system is denoted 
by {q;p}, where q = {q lx , q iy , q lz , ■ ■ ■ , qNx, qNy, qNz} 
and p = {pix,Piy,Piz, ■ ■ ■ ,PNx,PNy,PNz} represent the 
configuration and momentum spaces, respectively. The 
phase space contraction factor, A = + f^, depends on 
the detail dynamics of the system. 

Assume that the classic system contacts with a ther- 
mal bath at constant temperature T, and that it stays 
at an equilibrium state for time t < t\. Take a micro- 
scopic state A\ corresponding to this equilibrium state. 
From time t\ to t' 2l we switch on an external field denoted 
by a parameter A varying from Ai to A2, and drive the 
system out of equilibrium. After a sufficient relaxation 
with fixed A2 , the system arrives at the other equilibrium 
state at time t%. Correspondingly, the microscopic state 
evolves to A 2 . From time t\ to t 2l the entropy produc- 
tion function along the microscopic path <y(t) linking the 
states A\ and A 2 is expressed as [jj 

S [ 7 (*)] = Hh/h) - [' A[ 7 (t)]A, (1) 

where j\ and /2 are the equilibrium distribution func- 
tions at time t\ and £2, respectively. One can prove the 
FT, Pf(s)/pr(— s) — e s , wherepF(s) andp_R(s) represent 
the probability distributions of the entropy production 
function taking value s along the microscopic path "/(t) 
and its time-reversal path, respectively. If averaging e~ B 
for all paths beginning from all microscopic states cor- 
responding to the macroscopic equilibrium state at time 
t\, we have 

<e- s ) = J e- s p F (s)ds = J p R (-s)ds = 1. (2) 

This is nothing but the Kawasaki identity or Hatano-Sasa 

equality [2ll27|. 

If taking canonical distributions for the initial state 

at time t\ and the final state at time t%, we have /1 = 
e 0{F x -H x ) and j 2 = e P(F 2 -H 2 )^ where Fi and p 2 are the 

free energies of the system at time t\ and t 2 while H\ and 
H 2 are the Hamiltonians of the system at time t\ and t 2 . 
Assume that the effective dynamics of the system can be 



expressed as [23j 

q„ = dH/dp n , (3) 
p„ = -dH/dq n - a[7(t)]p„, (4) 

where a[7(t)] is the thermostat multiplier |2^| ensuring 
the kinetic temperature of the system to be fixed at T, 
and it reflects the heat exchange between the system and 
the thermal bath. H is the A-dependent Hamiltonian. 
Under the above dynamics, the phase space contraction 
factor is derived as A[7(i)] = — 3Na[y(t)] and its integral 
from time t\ to ti is just the entropy production induced 
by the heat (g[7(£)]) absorbed by the system from the 
thermal bath along the microscopic path j(t) linking the 
states A\ and A 2 , i.e., 




A[<y(t)}dt = 0q[ 7 (t)}. (5) 



This equation is crucial to the derivation of JE from the 
FT. Thus the entropy production function, Eq. JQ), is 
transformed into 

s[ 1 (t)]=0(wh(t)}-AF), (6) 

where w[7(i)] = H 2 — H 1 — g[7(t)] (microscopic energy 
conservation) is the work performed on the system along 
the microscopic path 7(i). AF = F 2 — F\ is the free 
energy change of the system from time t\ to t 2 . Assume 
that u>[7(£)] takes value w when s[7(t)] has value s, and 
notice that there is no work from time t' 2 to t 2 because 
the parameter A is unchanged at this time interval. From 
Eqs. (01 and JBJ) we easily arrive at the JE, 

( e -*») = e ^ AF . (7) 

We emphasize again that Eq. (JSJ), the time integral of 
the phase space contraction factor exactly expressed as 
the entropy production resulting from the heat absorbed 
by the system from the thermal bath, is the crucial point 
in the derivation of the JE from the FT. Remember that 
the phase space contraction factor depends on the mi- 
croscopic dynamics. If the dynamics satisfies the con- 
dition of adiabatic incompressibility |23| . i.e., the phase 
space contraction factor depends merely on the thermo- 
stat multiplier, Eq. © holds and the JE is a natural 
corollary of the FT. 

III. GENERALIZED JE 

Now, we check whether the condition of adiabatic in- 
compressibility holds for the systems mentioned in the 
above two experiments. 

For the first experiment, because the initial velocity 
distribution obeys the Maxwell distribution, some parti- 
cles with velocity larger than v p will strike the piston and 
then bounce fully but v p is unchanged because the mass 
of the piston is much larger than the total mass of the 
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particles that collide with it. In each bounces, the piston 
will do a small work on the gas system. Lua et al. have 
proved that the mean work is vanishing but the JE still 
holds for v p — > 
expressed as [25 



|l8j. The effective dynamics can be 



q« = Pn/™ + R--q«, 
p„ = -R • p n — a[j(t)]p n , 
V = Vr, 



(8) 
(9) 
(10) 



where m is the mass of each particle, the matrix R = 
((f, 0,0)', (0,0,0)', (0,0,0)*), and f the time-dependent 
volume expansion ratio with vanishing value except at 
the time interval between t% and t' 2 . Here the volume V 
plays the role of the parameter A in Evans' derivation of 
JE from the FT. The phase space contraction factor is 
found to be A [7(4)] = -3iVa[ 7 (i)]. Hence Eq. J5J) as well 
as the condition of adiabatic incompressibility still holds 
and so the JE is valid in this experiment. 

The effective term — R • p„ in Eq. J^J reflects the col- 
lisions between the piston and particles. For the second 
experiment, the volume expansion has no direct effect on 
the momentum of the particles. Thus the term — R • p„ 
should be removed, and Eq. should be replaced by 



Pri = -a[7(*)]Pr, 



(11) 



but Eqs. (JSJ and (|10|) are kept intact in the effective dy- 
namics of the second experiment. Consequently, we ob- 
tain the phase space contraction factor A [7(f)] = Nr — 
3Na[y(t)] and its integral from time t\ to 



A[ 7 (t)]dt = Nhx{V 2 /V x ) + 0qfr(t)], 



(12) 



ti 



where we have used J* 2 fdt — Jy^ d(lnV) = ln(V2/Fi) 

and - J t * 2 3Na[j{t)]dt = Pq[y(t)}. Obviously, Eq. © as 
well as the condition of adiabatic incompressibility does 
not hold in this experiment and so the JE fails. How- 
ever, following the derivation from Eq. 10 to Eq. J7J), 
and replacing Eq. (J5J by Eq. i|12[) . we obtain a general- 
ized equality beyond JE: 



-0[w-NT\n(V 2 /Vi)]\ 



-PAF 



(13) 



Because w = and AF = -ATTln(F 2 /Fi) in the sec- 
ond experiment, the above equation holds although the 
original form of JE fails. 

Through the above discussions, we know that the JE 
holds in the first experiment but fails in the second one. 
The underlying reason is that these two experiments have 
different microscopic dynamics: One satisfies the condi- 
tion of adiabatic incompressibility but another does not. 
Especially, the second experiment suggests that a more 
general form of JE should exist, and we would like to 
consider this possibility. Enlightened by Eq. (fL2"|> . we 
divide the time integral of the phase space contraction 
factor into two parts: One is /3q[y(t)], the entropy pro- 
duction resulting from the heat absorbed by the system 



from the thermal bath; Another is the entropy induced 
by the change of the external parameter and expressed 
as (3a. That is, 



A[j(t)]dt = (3(q[l(t)} + t). 



(14) 



Following the derivation from Eq. (JSJ to Eq. (J7J and re- 
placing Eq. (J3J) by Eq. (|14|) . we arrive at a generalized 
JE: 



-/3(w-<t) 



which is transformed into 



-0w\ 



-PAF 



-/3(AF+<r) 



(15) 



(16) 



if a depends only on the value of the external parame- 
ter at time t\ and i2, but not explicitly on the micro- 
scopic pathes. We conjecture that a = for most macro- 
scopic system and then Eq. I|15|) is degenerated into the 
original JE. a 7^ only in some very special systems 
and correspondingly the original JE fails. For example, 
fj = NT \n(V2/Vi) in the second experiment. 

Noticing that the mass of the piston in the first exper- 
iment is infinitely large. If it has an infinitesimal value, 
this system is equivalent to that in the second experiment 
because the collisions between particles and the piston 
has no effect on the momentum of the particles such that 
the particles do not feel the existence of the piston. It 
is interesting to discuss the intermediate case between 
the above two limits. Let us consider the third gedanken 
experiment where the experimental setup is the same as 
the first one except the mass M of the piston is finite. At 
time t\, we remove the pins Pi and P 2 , and the gas will 
push the piston to the right wall of the container. Once 
the piston contacts with the wall, it adheres to the wall 
without bounce. After a long time relaxation, the system 
arrives at an equilibrium state at time t% . When we write 
the effective dynamics, Eq. should be replaced by 



.gR- p n - a[7(t)]p„, 



(17) 



but Eqs. JSJl and (|rU|) are unchanged, where g is a function 
of m and M taking values between and 1. g may also 
depend on N and V2/V1 because the equations of motion 
are just the effective ones. With this dynamics, we obtain 
a = (1 - g)NT\n(V 2 /V 1 ) from Eq. |Q3J}. Thus Eq. JT5J 



gives 



ln( 



-0W\ _ 



)=gNMV 2 /V 1 ). 



(18) 



In order to recover the former two experiments, g must 
satisfy g — » 1 for M — > 00 and g — > for M — * 0. To 
determine g, we do numerical simulations for ideal gas 
with different N (from 1000 to 10000), M/m (from 0.2 
to 1000), and V 2 /Vi (from 1.1 to 1.9), and calculate g 
by Eq. (|18fl . To obtain the ensemble average (e~^ w ), we 
take 500 systems j2^| with different initial microstates 
corresponding to the same macroscopic equilibrium state. 
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We find that g depends only on the combined variable 
x = M/[mN\n(V2/Vi)]. The relation between g and 
x is shown in Fig. For very small x, the numerical 
data (the inset of Fig. |2J) can be fit well by a line \ng = 
0.93 + 0.88 In x. That is, g has the asymptotic form g ~ 
(2.87x) 88 for x -> (corresponding to M -> 0). Based 
on this asymptotic form, noting that <7 — ► 1 for a; — ■> oo 
(corresponding to M — > oo), we conjecture that g has the 
form 



(2.87a;) 1 ' 
1 + (2.87x)" 



0.88/v 



(19) 



Our numerical data is indeed fitted well by this form. The 
fitting curve is the dash line in Fig. [21 with the parameter 
v = 0.53. We use this fitting parameter and Eq. i|19|) 
to predict g = 0.8795 for M/m = 4000, N = 1000 and 
V2/V1 = 1.1, which is quite close to the value 0.8846 ob- 
tained from the numerical simulations. This fact implies 
that our conjecture is reasonable although we cannot in- 
tuitively figure out the physical meaning of the numbers 
2.87, 0.88 and v = 0.53 in Eq. fTS} . 



push Ceo to the new position, and one can measure the 
velocity of Cgo when it arrives at the new position and 
calculate the corresponding work. Repeat this process 
for many times and calculate the value of — ln(e~P w ) / /3. 
Comparing this value with the free energy obtained from 
theoretical calculation, one can obtain the value of a. If 
a 7^ 0, the JE is violated and a generalized JE should 
exist. 

Long (1 0,1 0)SWNT 




FIG. 3: Schematic figure of the experimental setup (in vac- 
uum). 



IV. DISCUSSION AND CONCLUSION 
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FIG. 2: Numerical results and fitting curve for the relation 
between g and x where x represents M/[mN ln(V2/Vi)]. The 
squares come from numerical simulations. The result for small 
x < 0.01 is magnified in the inset of the figure. 



For the macroscopic gas system except for the case 
in the second experiment, we have in general M 3> 
mN ln(V 2 /Vi), so g ~ 1. Hence a = and Eq. {T^J 
is degenerated into the JE. Therefore the departure from 
the JE should occur at the small scale system with fi- 
nite M but still large enough N. For example, our result 
might be verified for the inert gas in a very long single- 
walled carbon nanotube (SWNT) as shown in Fig.[3J One 
end of the nanotube is closed while another is opened, 
and a buckyball Cgo is put in it as a piston. Select the 
proper nanotube, for example (10,10) nanotube, and the 
gas with large radius, for example Ar, such that Cgo can 
prevent the gas from escaping from the interstice between 
Ceo and the nanotube. A small SWNT can be used to fix 
the initial position of Ceo- At some time, pull the small 
SWNT outward to another position quickly. The gas will 



It is useful to discuss some questions before concluding 
this paper. 

(i) If the volume of the thermal system is fixed but 
other parameter varies such as in the single molecule 
mechanical experiments 0, JE always holds be- 
cause the effective dynamics can be expressed as Eqs.© 
and l|5|). Generally speaking, one may not construct a 
parameter-dependent Hamiltonian and express the effec- 
tive dynamics as Eqs.® and |J2J| if the volume changes. 
JE may be violated in this case. However, JE still holds 
if one can control the ratio of volume change because the 
the particles fully bounce when they collide with the pis- 
ton in this case, i.e., controlling the ratio is equivalent to 
M -> 00. 

(ii) The derivation of JE from FT does not require the 
thermostat for the whole process from t\ to t%. For t < t\, 
the system is at an equilibrium state in contact with a 
thermal bath at temperature T. The external field is 
switched on from time t% to t' 2 . The contact with a ther- 
mal bath is unnecessary in this stage (i.e. no heat ex- 
change, q = 0). After that, the external field is fixed and 
let the system contact with the same thermal bath. At 
time £2 j the system reaches the equilibrium state through 
a long time relaxation. The thermostat is merely required 
in this stage (from t' 2 to t 2 ). In fact, this requirement is 
the same as the original proof of JE [6j , which ensures us 
to calculate (e~^ w ) easily from numerical simulations. 

(iii) Equations P)l-l(TUj). l|H)l. and ljTT|) are the effective 
dynamics of thermal system with the volume changes. 
Here the "effective" means that the dynamics is not the 
real microscopic motion (of course, the real motion for 
each particle still abides by Newtonian laws), while it is 
the image mapping from the real motion and can give the 
correct thermodynamic properties of the system through 
Molecule Dynamics Simulations [2j|. The effective term 



5 



— R • p n in these equations reflects the collisions between 
the piston and particles. If the particle fully bounces 
(i.e. M — ► oo), the coefficient before this term is 1. If 
the particle do not bounce (i.e. M — > 0), the coefficient 
is 0. For finite M, the coefficient should intermediate 
between and 1. We use g to express this coefficient in 
Ea. (|17|l . In the second experiment, V = V\ for t < t\, 
but V = Vi for t — t' 2 — ■> t~l . Thus f is a 5- function, which 
implies q„ changes discontinuously. This is impossible for 
real dynamics but permitted in the effective dynamics. 
Our numerical simulation is performed for real dynamics 
(Newtonian mechanics) and the results reveal that g — » 
1, a — ► thus JE holds for M — > oo, and that g — » => 
ct = iVrin(F 2 /Fi) thus Eq. JEJl holds for M -> 0. That 
is, the numerical results obtained from the real dynamics 
are the same as those derived from the effective dynamics 
Eqs. JSJ (|10|l and (|ll|l . which suggests that the effective 
dynamics is consistent with the real dynamics and our 
argument that the term — R • p„ reflects the collisions 
between the piston and particles is reasonable. 

In summary, we have pointed out that the crucial point 
in the derivation of the JE from the FT is that the time 
integral of the phase space contraction factor is exactly 
expressed as the entropy production resulting from the 



heat absorbed by the system from the thermal bath, i.e. 
the dynamics of the system satisfies the condition of adi- 
abatic incompressibility. For the system violating this 
condition, a more general version of JE, Eq. (|15|l . exists. 
In the future, it is interesting to find some real systems 
which makes a ^ 0. Deriving the analytic expression 
of the quantity a is another challenge because a might 
be system-dependent. These researches will enhance our 
understanding to non-equilibrium statistics. 
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